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ABSTRACT 


Previous  investigations  have  demonstrated  the  importance  of  the  effect 
of  linearly  varying  axial  or  in-plane  loading  on  the  vibration  characteristics 
of  beams  and  flat  plates.  It  has  already  been  established  that  the  problem 
reduces  to  solving  for  the  eigenvalues  of  a  fourth  order,  variable  coefficient 
differential  equation  that  can  not  be  solved  in  closed  form.  Beginning  with  a 
variational  representation  of  the  eigenvalue  problem,  methods  are  discussed 
by  which  both  upper  and  lower  bounds  for  the  eigenvalues  may  be  formed.  The 
true  eigenvalues  may  thus  be  estimated  as  being  bracketed  by  the  upper  and 
lower  bounds  which  are  shown  to  approach  each  other.  The  bounds  for  the 
eigenvalues  may  also  be  estimated  by  an  averaging  procedure  which  may  or  may 
not  compare  favorably  with  the  true  values  depending  on  the  values  of  the  load 
ing  parameters.  Finally,  numerical  values  for  upper  bounds,  lower  bounds,  and 
average  lumped  end-load  eigenvalues  are  computed  on  an  IBM  7090  Computer. 
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NOMENCLATURE 


Differential  operator  of  loaded  beam 

Eigenvectors 

Constants 

Modulus  of  Elasticity 

Natural  frequency  of  vibrating  beam 

Moment  of  inertia 

Class  of  admissible  functions  in  elastic  stability  problems 

Class  of  admissible  functions  in  vibration  problems 

Length  of  the  beam 

Constant  end  loads 

Function 

Function 

Axial  coordinate 

Distributed  axial  load  parameter 

Critical  axial  load 

Ratio  of  end  load  to  total  distributed  load 

Separation  constant 

Function 

Eigenvalue 

Upper  bound 

Nondimens ional  axial  variable 

Density  per  unit  of  length 

Mode  shape,  dependent  deflection  variable 


I.  INTRODUCTION 


In  recent  years  much  attention  has  been  given  to  the  effect  of  lin¬ 
early  varying  axial  or  in-plane  loads  on  the  vibrational  characteristics  of 
beams  and  plates.  This  topic  is  of  particular  interest  in  aerospace  appli¬ 
cations  where  inertia  and  friction  drag  forces  manifest  themselves  as  axial 
or  in-plane  loads,  A  detailed  formulation  of  the  problem  is  the  subject  of 
a  prior  NASA  report  by  authors  and  is  the  subject  of  considerable  liter¬ 
ature  (see  ref.  2  through  16). 


Formulation  of  the  Problem 

As  described  in  references  (1)  and  (2),  the  eigenvalue  problem  for 
both  the  beam  and  the  rectangular  plate  may  be  resolved,  under  certain  re¬ 
strictions,  to  a  solution  of  the  ordinary  differential  equation 


(1) 


and  the  boundary  conditions 
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In  view  of  the  definition  of  the  parameter  ^  ,  it  is  clear  that  for  a 
given  compressive  distributed  load  (jj  3  the  follov/ing  cases  may  occur: 


1)  0  *  beam  is  entirely  in  compression 

beam  is  partly  in  tension  and  partly  in  compression 

"1^/1  ’  entirely  in  tension  since  the  tensile  and 

load  Pj  is  larger  than  the  total  distributed  load  L. 

In  the  last  case,  the  problem  of  elastic  stability  does  not  exist* 

The  determination  of  mode  shapes  and  natural  frequencies  involves  the 
solution  of  the  differential  eigenvalue  problem  defined  by  eqs,  (1)  and  (2). 
Variational  techniques  (  )  (  )  finally  resolve  this  to  obtaining  solutions  to 
the  variational  principle 


(4) 


where  K  is  the  class  of  functions  constituting  the  domain  of  definition  of 
the  operator  A,  and,  hence,  satisfying  both  the  prescribed  and  the  natural 
boundary  conditions,  and  <  >  denotes  the  inner  product  between  tt70 

functions  u,v  ,  where 


and 


(5) 


Equation  (4)  may  be  characterized  by  Courant's  maximum-minimum  char¬ 
acterization  (ref.  18,  Chap.  Ill)  given  by 
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where  <j>  and  belong  to  K  ,  where  is  the  class  of  ad 
missible  functions  required  to  satisfy  only  the  prescribed  boundar 


y  only  the  prescribed  boundary  conditions 


*  This  functional  is  known  as  Rayleigh's  quotient. 
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In  resume,  the  situation  is  as  follows:  if  p  is  such  that  buckling 
may  occur,  there  exists  for  the  given  falue  of  p  a  critical  value  of  this 
distributed  axial  load  parameter,  ,  for  which  the  beam  is  unstable  and 

the  potential  energy  is  equal  to  zero.  For  any  value  of  cv  less  than  , 

the  potential  energy  is  positive,  and  the  beam  has  discrete  natural  frequen¬ 
cies  whose  square  are  proportional  to  the  eigenvalues  of  the  operator  A, 

where 


These  eigenvalues  are  assumed  to  be  ordered  in  the  non-decreasing 
sequences 

. 


The  eigenfunctions  corresponding  to  distinct  eigenvalues  are  mutually  orthog 
onal,  and  correspond  to  the  mode  shapes  of  the  beam.  For  a  given  value  of 
p  ,  as  a  increases,  the  numerator  of  the  Rayleigh  quotient  decreases  and 
the  eigenvalues  decrease.  Buckling  occurs  v/hen  or  becomes  equal  to  , 

for  which  the  first  eigenvalue  goes  to  zero. 

In  the  next  section,  we  review  the  methods  used  in  this  work  to  ob¬ 
tain  approximate  solutions. 


II.  BOUNDS  FOR  EIGENVALUES 


There  appear  in  the  literature  many  m’ethods  for  finding  bounds  for 
eigenvalues.  Upper  bounds  are  usually  found  without  too  many  difficulties 
by  the  Rayleigh-Ritz  method.  Lower  bounds  present  considerably  more  diffi¬ 
culties,  and  it  can  be  said  that  no  method  having  the  generality,  simplicity 
and  success  of  the  Rayleigh-Ritz  method  exists  for  the  computation  of  lower 
bounds.  The  most  suitable  method  usually  depends  on  the  problem  at  hand. 

In  this  section,  we  review  briefly  the  methods  used  in  this 

work  in  the  calculations  of  approximations  to  eigenvalues.  They  are  the 
Rayleigh-Ritz  method,  the  method  of  Kato,  and  the  method  of  intermediate 
problems  of  Weinstein  and  Aronszajn,  with  some  modifications  introduced  by 
Bazley  and  Fox, 


A.  The  Rayleigh-Ritz  Method 


The  Rayleigh-Ritz  method  for  numerical  computations  of  approximations 
to  eigenvalues  has  been  used  extensively  and  with  great  success  in  the  liter¬ 
ature.  Consequently,  it  will  only  be  outlined  briefly  here. 

.The  basic  idea  of  the  method  consists  in  determining  the  stationary 
values  of  the  Rayleigh  quotient,  not  over  all  admissible  functions  u,  but 
only  over  the  linear  manifold  spanned  by  an  arbitrary  set  of  n  linearly  in¬ 
dependent  functions  satisfying  the  boundary  conditions  of  the 

operator  A.  The  problem  then  consists  in  finding  the  functions  u  of  the  form 


JUlx 


A 


(8) 


finding  the  constants  Cj|_,  making  the  Rayleigh  quotient  stationary, 
and  the  stationary  value  of  the  quotient.  Substitution  of  Equation (8) into 
Rayleigh's  quotient  yields 
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which  is  the  ratio  of  two  quadratic  forms  in  the  n  real  variables  C-i ,  Co, 
...Cn.  Its  stationary  values  can  be  obtained  by  finding,  for  instance 
the  stationary  values  of  the  quadratic  form  in  the  numerator,  subject  to 
the  auxiliary  condition  that  the  denominator  be  equal  to  one,  and  using 
the  method  of  the  Lagrange  undetermined  multiplier.  The  result  is  the 
general  matrix  eigenvalue  problem 


[<  AU,  Att-pJ  [cj  J  A  (10) 


See,  for  instance,  references  17,  18,  and  19. 

The  functions  u^j^  are  often  called  coordinate  functions. 
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Since  the  class  of  admissible  functions  was  res^icted  to  the  pnite 
dimensiohal  manifold,  it  follows  that  the  eigenvalues  ^  are  upper-bounds 
for  the  eigenvalues  of  A,  i.e.> 


Furthermore,  it  follows  that  as  n  increases,  the  upper  bounds  are  improved, 
or  at  least,  not  worsened. 

From  a  computational  standpoint,  it  is  advantageous  to  choose  mut¬ 
ually  orthogonal  coordinate  functions  to  avoid  the  solution  of  a  genera 
eigenvalue  problem.  Also,  Equation(9)may  be  written  as  in  Equation (6 )with 
functions  |Ui\  required  to  satisfy  only  the 
tions.  This  point  is  discussed  in  detail  in  references  18  and  19.  The  co 
ordinate  functions  utilized  in  this  work  satisfy  both  the  prescribed  and  the 
natural  boundary  conditions,  as  will  be  seen  later. 


B,  The  Method  of  Kato 

The  Rayleigh-Ritz  method  described  above  furnishes  upper  bounds  for 
eigenvalues.  The  results,  particularly  for  the  first  eigenvalue,  are  usually 
in  agreement  with  the  exact  eigenvalues  for  the  cases  where  the  latteicn 
be  obtained.  However,  in  general,  the  question  regarding  the  closeness  of 
these  bounds  to  the  true  values  remains  unanswered,  although  in  some  instances, 
an  estimate  of  the  error  is  possible.  One  way  of  determining  how  good  the 
approximations  are  is  to  compute  also  lower  bounds.  If  ^hese  turn  out  close 
to  the  upper  bounds,  the  question  is  essentially  answered.  The  method  of 
j^3t.o(22)^  which  is  an  extension  of  Temple's  method,  furnishes  lower  bounds, 
provided’ rough  estimates  of  the  sought  eigenvalues  are  known.  This  is  out¬ 
lined  by  Freidman  P*  ^1^), 


C*  The  Method  of  Intermediate  Problems 

The  methods  described  in  the  preceding  t\Jo  sections  furnish*  upper  and 
lower  bounds  for  eigenvalues.  In  both  methods,  the  quality  of  the  results 


See,  for  instance,  reference  21,  p.  336. 
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depends  strongly  on  how  well  the  trial  functions  approximate  the  eigenvectors 
of  the  operator.  Hence,  both  methods  may  require  considerable  ingenuity  in 
the  selection  of  the  trial  functions.  Furthermore,  for  different  sets  of 
trial  functions,  there  is  little  prior  knowledge  of  which  set  will  give  the 
best  results.  For  these  reasons,  it  is  in  order  to  consider  also  another 
method  for  the  computation  of  the  lower  bounds.  The  method  used  here  is  the 

method  of  intermediate  problem,  which  presents  the  advantage  that  the  bounds 
can  be  improved. 

Quite  a  few  years  back,  Weinstein  introduced  the  method  of  in¬ 

termediate  problems,  which  gives  improvable  lower  bounds  by  changing  the 
boundary  conditions  of  differential  operators.  Briefly,  the  method  consists 
in  relaxing  the  boundary  conditions  to  obtain  a  solvable  problem,  the  base 
problem,  whose  eigenvalues  give  rough  lower  bounds  for  the  eigenvalues  of 
the  given  problem.  A  sequence  of  intermediate  problems  linking  the  base  pro¬ 
blem  to  the  given  problem  is  then  introduced.  These  are  such  that  they  can 
be  solved  In  terms  of  the  base  problem,  and  that  they  give  improved  lower 
bounds.  The  details  of  the  procedure  are  exposed  in  references  17  and  25. 

^  1951,  Aronszajn  ''pointed  out  that  a  base  problem  can  be  ob¬ 

tained  by  changing  the  differential  operator,  and  Indicated  the  method  of 
construction  of  the  intermediate  problems.  The  solution  of  these  intermed¬ 
iate  problems  requires  the  determination  of  the  poles  and  the  zeroes  of  a 
meromorphic  function  given  in  its  partial  fractions  representations.  From 
a  c^putational  standpoint,  the  determination  of  the  zeroes  present  many 
difficulties  which  have  been  removed  in  a  dissertation  by  Bazley  (^7) 

Barley  and  Fox  (28-33)^  authors  have 

applied  their  method  to  the  determination  of  the  eigenvalues  of  Schrodlnger's 
equation  and  Mathieu's  equation  with  excellent  results. 

A  more  detailed  resume  of  the  Method  of  Rato  and  the  Method  of  In- 
te^ediate  Proble^  is  given  in  Reference  (2).  Reference  (2)  also  describes 
specific  application  to  the  simply  supported  beam  and  the  beam  with  built- 
in  ends.  These  procedures  are  not  particularly  difficult  in  principle,  but 
the  calculations  Involved  are  somewhat  laborious 


D.  Lumped  Constant  End  Load  Approximation 

1  j  ^PP’^o^^^^tion  to  the  response  of  beams  with  distributed  axial 

load  may  be  accomplished  by  replacing  the  distributed  load  and  its  reaction 
with  equal  and  opposite  average  end  loads.  This  results  in  an  ordinary  lin- 
ear  differential  equation  with  constant  coefficient  which  may  be  solved  ex¬ 
actly  in  terms  of  trigonometric  functions.  A  comparison  of  the  eigenvalues 
calculated  in  this  manner  is  made  with  the  upper  and  lower  bounds  in  the 
section  on  Results  and  Discussion. 
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Ill,  RESULTS  AND  DISCUSSION 


Following  the  methods  described  above,  upper  and  lower  bounds  for  the 
eigenvalues  of  the  simply  supported  and  clamped  beam  were  calculated  on  an  IBM 
7090  Computer  in  the  Computation  and  Data  Processing  Center  of  the  University 
of  Pittsburgh,  The  results  are  displayed  in  Tables  I  and  II  and  Figures  1,  2, 
and  3,  Uppar  bounds,  lower  bounds  and  lumped  end- load  eigenvalues  are  dis¬ 
played  for  a  wide  range  of  loading  parameters  a  and  p  . 


A,  Simply  Supported  Beam 

The  bounds  for  the  first  five  eigenvalues  of  the  simply  supported  beam 
are  presented  in  Table  I,  To  facilitate  the  comparison  between  the  Rayleigh- 
Ritz  upper  bounds  and  the  lower  bounds  by  the  method  of  intermediate  problems, 
the  ratio  of  their  difference  to  their  average  has  been  computed  and  is  also 
presented  in  Table  I.  Since  the  eigenvalues  of  a  simply  supported  beam  are 
easy  to  obtain,  it  is  interesting  to  compare  the  upper  and  lower  bounds  of  the 
eigenvalues  obtained  by  lumping  half  of  the  total  distributed  load  as  a  con¬ 
stant  load  at  each  end.  These  results  are  also  included  in  Table  I, 

Analysis  of  the  results  in  Table  I  indicates  that  the  Rayleigh-Ritz 
upper  bounds  and  the  lower  bounds  by  the  method  of  intermediate  problems  re¬ 
main  close  over  the  whole  range  of  axial  loadings.  This  is  particularly  true 
for  the  first  eigenvalue.  Only  when  the  beam  is  extremely  close  to  buckling 
does  the  relative  error  increase  greatly  as  a  result  of  the  smallness  of  the 
eigenvalues.  For  eigenvalues  of  order  higher  than  one,  the  error  is  slightly 
higher,  but,  if  necessary,  it  could  be  reduced  by  considering  higher  intermed¬ 
iate  problems. 

The  lower  bounds  for  the  first  eigenvalue  by  the  method  of  Kato  re¬ 
main  close  to  the  upper  bounds  for  moderate  loading,  but  drop  off  considerably 
at  the  loading  increases.  Perhaps,  this  effect  might  be  attributed  to  the  fact 
that  as  the  first  eigenvalue  approaches  zero,  the  choice  arbitrary  trial  varia¬ 
tions  becomes  more  and  more  critical.  For  higher  eigenvalues,  this  selection 
is  not  as  critical,  and  consequently,  the  lower  bounds  remain  close  to  the 
upper  bounds.  However,  in  the  cases  where  the  beam  can  not  become  elastically 
unstable,  the  Kato  lower  bounds  eventually  decrease  as  the  loading  becomes  very 
large,  and  no  explanation  for  this  behavior  can  be  offered. 

The  eigenvalues  of  the  beam  with  lumped  constant  end  load  are  remark¬ 
ably  close  to  those  of  the  beam  with  distributed  load  for  compressive  end 
thrusts,  i,e,,  forp>o  .  Ebr  negative  p  ,  the  results  are  quite  far  apart. 

In  particular,  forp=-5,  the  beam  with  distributed  axial  load  may  become 
elastically  unstable,  while  the  beam  with  lumped  load  can  not  buckle,  because 
its  net  thrust  is  zero.  Consequently,  extreme  care  should  be  exercised  in  the 
lumping  of  the  loads  when  they  are  of  opposite  signs. 
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The  effect  of  the  axial  loads  on  the  first  frequency  of  the  simply  sup- 
pori.ed  beam  is  shown  in  Figures  1  and  2.  Figure  1  represents  the  ratio  of  the 
first  frequency  of  the  loaded  beam  to  that  of  the  unloaded  beam  as  a  function  of 
tne  distributed  load  parameter  a  ,  as  obtained  by  Kato's  method  and  the  Ray- 
leigh-Ritz  method.  The  lower  bounds  of  the  method  of  intermediate  problems  are 
not  shown  because  their  curve  practically  coincides  with  the  Rayleigh-Ritz  curve 
.or  tne  scale  used.  The  curves  correspond  to  p =  o  .  Figure  2  also  represents 
the  ratio  or  the  fundamental  frequency  of  the  loaded  beam  to  that  of  the  unloaded 
beam  as  a  function  or  a  for  various  values  of  p  .  The  curves  were  obtained 

y  using  the  average  of  the  upper  bounds  and  lower  bounds  by  the  method  of  inter¬ 
mediate  problems. 


The  values  of  the  critical  axial  load 


a 


^  ^  .  -  -  c  given  at  the  intersection 

Oi.  ct-ie  rrequency  ratio  curve  with  the  horizontal  axis.  The  buckling  loads  ob¬ 
tained  from  graphs  having  a  larger  scale  thaii  that  of  Figure  2  compare  favorably 
with  the  exact  results  of  Tyler  and  Rouleau  po,-  p^O  ,  the  graphs  indi¬ 

cate  that0<5  18.7  Ei/l  v/hile  Tyler  and  Rouleau's  result  is  ot,  -  18  7filEl/L^ 
Fetp^l  O  ,  we  obtain  6.5  EI/L^  ' 

xvhile  the  exact  answer  is  =  6.519E.I/L^  and  for  p=  -.50  we  haveo<,»i  83  EI/L' 
against  the  exact  result  ol  82.8819  E I/L^  .  The  approximate  values  are  "certain¬ 
ly  close  enough  lor  engineering  application. 


B.  Clamped  Beam 

^  The  bounds  for  the  first  four  eigenvalues  of  the  clamped  beam  are  pre- 
sentea  in  Table  2.  The  ratio  of  the  difference  between  the  upper  bound  and  the 
corresponding  lower  bound  by  the  method  of  intermediate  problems  to  their  aver' 
age  has  also  been  computed.  The  eigenvalues  of  the  clamped  beam  carrying  a  con 
SLanL  enu  load  equal  to  half  the  total  distributed  load  and  the  constant  end 
loaci  are  also  presented  in  Table  II  to  indicate  for  what  values  of  the  loading 
parameters  this  lumping  is  acceptable.  ° 

Examination  of  the  results  indicate  the  following: 

i)  The  lower  bounds  by  the  method  of  intermediate  problems  are 
very  close  to  the  Rayleigh-Ritz  upper  bounds  for  all  eigen¬ 
values  and  for  the  v/hole  range  of  the  loading  parameters. 

ii)  The  lower  bounds  by  the  method  of  Kato  present  the  same  fea¬ 
tures  demonstrated  in  the  simply  supported  beam  calculations* 
whenever  the  loading  is  small,  the  bounds  are  fairly  good 
but  become  vjorse  as  the  loads  increase. 

iii)  The  eigenvalues  of  the  beam  with  lumped  end  load  are  fairly 
close  to  the  upper  bounds  for  moderate  loading,  particularly 
forp>0  .  For  negative  values  of  p  ,  they  can  be  quite 
remote  from  the  upper  bounds,  particularly  forP  for 

which  the  beam  with  distributed  axial  load  may  buckle  while 
the  beam  V7ith  lumped  end  load  can  not. 
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The  effect  of  the  axial  loads  on  the  first  frequency  of  the  clamped 
beam  are  shown  in  Figure  3,  which  represents  the  ratio  of  the  first  fre¬ 
quency  of  the  loaded  beam  to  that  of  the  unloaded  beam  as  a  function  of  the 
axial  load  parameter  a  for  various  values  of  p  , 


IV.  CONCLUDING  REMARKS 


Bounds  for  the  eigenvalues  of  a  simply  supported  and  a  clamped  beam 
carrying  linearly  distributed  axial  loads  have  been  presented.  The  main  dif¬ 
ficulty  in  problems  of  this  nature  arises  from  the  fact  that  the  governing 
differential  equation  has  a  varying  coefficient  which  usually  prevents  one 
from  obtaining  exact  solutions.  Upper  bounds  were  easily  obtained  by  the 
Rayleigh-Ritz  method.  Lower  bounds  by  the  method  of  Kato  were  also  easy  to 
obtain.  In  both  methods,  the  closeness  of  the  results  to  the  true  eigen¬ 
values  depends  on  the  quality  of  the  coordinate  functions.  It  appears  that 
for  moderate  loading,  the  eigenfunctions  of  the  unloaded  beams  were  good  co¬ 
ordinate  functions,  as  our  results  indicate. 

The  lower  bounds  computed  by  the  method  of  Intermediate  problems  were 
very  close  to  the  upper  bounds,  both  for  the  simply  supported  and  the  clamped 
beam.  The  modifications  introduced  by  Bazley  and  Fox  eliminate  the  computa¬ 
tional  difficulties  which  prevented  extensive  use  of  the  method  of  intermed¬ 
iate  problems . 

For  engineering  applications,  it  appears  that  lumping  the  axial  loads 
gives  eigenvalues  that  are  larger  than  the  true  eigenvalues,  and  that  care 
must  be  exercixed  whenever  the  distributed  load  and  the  constant  end  thrust 
are  of  opposite  signs.  In  this  case,  the  buckling  loads  predicted  by  the 
lumped  end  load  problem  can  be  quite  remote  from  the  actual  critical  loads. 

The  present  research  could  be  extended  to  the  consideration  of  beams 
with  other  boundary  conditions,  closer  determinations  of  the  buckling  loads, 
and  the  methods  used  here  can  be  applied  to  other  problems  giving  rise  to 
differential  equations  with  variable  coefficients,  such  as  in  the  problems 
of  the  determination  of  natural  frequencies  and  buckling  loads  of  beams  of 
varying  cross  sections,  plates  with  varying  in-plane  loads,  and  plates  of 
non-uniform  thickness,  to  mention  a  few.  Information  of  this  nature  would  be 
valuable  to  designers,  particularly  in  the  Aerospace  industry. 
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Effect  of  Axial  Loads  on  the  Fundamental  Frequency  of  a  Simply  Supported  Beam 
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4  34  574.8  34  344.0  34  537.7  0.107  34  579. 
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